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Current neutrino oscillation data indicate the existence of two large lepton mixing angles, while 
Kobayashi-Maskawa matrix elements are all small. Here we show how supersymmetric SO(IO) with 
extra chiral singlets can easily reconcile large lepton mixing angles with small quark mixing angles 
within the framework of the successful Fritzsch ansatz. Moreover we show how this is fully consistent 
with the thermal leptogenesis scenario, avoiding the so-called gravitino problem. A sizeable asym- 
metry can be generated at scales as low as possible within the leptogenesis mechanism. We present 
our results in terms of the leptonic CP violation parameter that characterizes neutrino oscillations. 
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I. INTRODUCTION 

An endemic difficulty in unified gauge models is the understanding of fermion masses. Thanks to the 
brilliant series of experiments which led to the discovery of neutrino oscillations over the past ten years or 
so, we now have an improved knowledge of neutrino masses and mixings [ij. The maximal atmospheric 
mixing angle and the large solar mixing angle both came as a surprise, at odds with naive unification- 
based expectations that lepton and quark mixings are similar. Indeed, although not mandatory, within 
unification models a similar structure for quark and lepton mixing is far easier to account for than what 
is currently established by observation. Several approaches have been considered in the literature to 



iQ,flaa. 



circumvent this problem, involving various types of extensions of the multiplet content 

A virtue of seesaw models is that they bring the possibility of decoupling the angles in the quark and 
lepton sectors by having additional Majorana-type terms for the latter Moreover, seesaw models may 
account for the observed cosmological baryon excess through the elegant mechanism called leptogene- 
sis However in minimal supergravity models, with 7713/2 ^ 100 GeV to 10 TeV the supersymmetric 
type-I seesaw scenario leads to an overproduction of gravitinos after inflation [i3,[ii| unless the reheat 
temperature is restricted to be much lower than that required for successful leptogenesis One pos- 
sible way out is resonant leptogenesis 13 1 as indeed suggested in Ref. Another alternative has been 



considered in Ref. [isl requires going beyond the minimal seesaw by adding a small R-parity violating 
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term in the superpotential. 

Here we consider an alternative approach to supersymmetric SO(IO) unification that involves mainly 
an extension of the lepton sector and a novel realization of the seesaw mechanism already discussed 



previously [161 . Il7l |. We show how the observed structure of lepton masses and mixings fits together, 
thanks to the presence of the extra states, with the requirements of a predictive pattern of quark mixing, 
encoded within the successful Fritzsch ansatz 18]. Moreover, we show explicitly how such extended 
supersymmetric SO(IO) seesaw scheme provides a realistic scenario for thermal leptogenesis, avoiding the 
so-called gravitino problem. Our approach is phenomenological in sipirit. Hence we do not seek to derive 



18, 



2C, 
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22, 



a full-fledged unified theory of flavour incorporating a specific symmetry, such as 
to the unified gauge group level. 

The paper is organized as follows. In Sec. HI] we discuss the limitations of minimal supersymmetric 
type-I seesaw in naturally providing a framework for thermal leptogenesis in agreement with bounds on 
the reheat temperature after inflation, in Sec. IIIII we recall the basic features of the considered model, 
while ansatzes for the coupling matrices are discussed in Sec. IIVI In Sec. |V]we discuss the calculation of 
the final baryon number asymmetry from the decay asymmetry of the lightest singlet. There we discuss 
how we can use results obtained previously on production and washout factors in the seesaw type-I case 
also in our extended seesaw model. We show that leptogenesis scenario is consistent both with the reheat 
temperature constraint as well as with fitting the lepton and quark mixing angles within the successful 
Fritzsch ansatz. 



II. MINIMAL SUPERSYMMETRIC SEESAW LEPTOGENESIS 

In the most general seesaw model the Dirac fermion mass terms arise from the couplings of the 16 
with the 10 and 126, while the diagonal entries in the seesaw neutrino mass matrix arise only from the 
couplings of the 16 with the 126. The mass matrix expressed in the basis is i^l, is given as j24j |. 



- / . u (i) 




where the 16 denotes each chiral matter generation while the 10 and 126 are Higgs-type chiral multiplets. 
Here Y, Yl and Yr denote the Yukawa couplings of the 10 and 126, respectively. These are all symmetric 
in flavour, the symmetry of Yl, Yr results from the Pauli principle, while that of Y follows from SO(IO). 
Minimization of the Higgs scalar potential leads to a vev seesaw relation 

(Ai) (A«) ^ , (2) 

consistent with the desired vev hierarchy (A/j) ^ ($) ^ (Al), where (<I>) is the standard vev and (Ai^ij) 
are the vacuum expectation values (vevs) giving rise to the Majorana terms. 

Before we display our proposed model it is useful to consider the case of the Minimal type-I seesaw [2^ . 
This consists in taking (A/,) — > in the above equations. Although not generally consistent with SO(IO) 
in the presence of the 126, this approximation is suitable for our purposes as the extension we will be 
discussing is tripletless, the 126 being replaced by 16's. 

The matrix M.^ is diagonalized by performing a a 6 x 6 unitary transformation U^, 

6 

^« = Y^{U,,).,ana U^MM^ = diag(m„ M,), (3) 
0=1 
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FIG. 1: Diagrams contributing to generate a lepton asymmetry of the universe. 



whose form is given exphcitly as a perturbation series in Ref. [26] . The effective mass matrix of the three 
light neutrinos is given as 

\2 



(Ar 



(4) 



The smallness of hght neutrino masses foUows dynamically from Eq. (j4|), since {Aji) is large. 

In order to reproduce the observed pattern of quark masses and CP violation we assume the complex 
matrix Y to have the form 



/ a 0\ 



Y 



c* d 



(5) 



with a, c complex and b, d real. Notice that the large observed lepton mixing angles can always be 
reconciled with the small quark mixing angles thanks to the presence of the coupling matrices 11,;?, which 
exist only for neutrinos (in the limit of unbroken D-parity one has Yl ~ Yu). However, in a unified theory 
where a flavour symmetry is assumed to predict the lepton mixing angles, it becomes a real challenge to 
account naturally also for the quark mixing angles. 

Now we turn to cosmology. One of the attractive features of seesaw models is that they open the 
possibility of accounting for the observed cosmological matter- antimatter asymmetry in the Universe 
through the leptogenesis mechanism Q . This requires the out-of-equilibrium decays of the heavy "right- 
handed" neutrinos Ni (see Fig. [T|) to take place before the electroweak phase transition, and the presence 
of CP violation in the lepton sector. The tree level and one loop diagrams for Ni decay that interfere in 
order to generate a lepton asymmetry of the universe are shown. 

These decays generate the asymmetry. 



r - 

^ NiL 



Fiv.L + r 



(6) 



NiL 



Here we recall the estimate of such asymmetry, as given in Ref. 27| . The relevant amplitudes follow from 
the initial Lagrangian 



C = X,jN,PlLjH + h.c.. 
Here it is sufficient just to sketch their structure. At the tree level we have 

while for the one-loop it is only necessary to evaluate the vertex contribution. 



(7) 
(8) 
(9) 
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where / is some complex loop function. Now tree level plus one loop will be proportional to 

\M%l + M],l\' ^ |A<VI'+>^^l>^]v*l+-M5V*l-M]vl (10) 
so that the asymmetry is 



Now using Eq. ([8]) and Eq. ([9]) we can write 



,2 
NLl 



(11) 



ocIm(/) ^ , (12) 



J2j Ki^k 

Finally we get, in the case of hierarchical right-handed neutrinos 



1 



X] 1"! (^Jl^ml'^mfe'^^i*:) 

(13) 

fe 

Note that in this limit the contribution of the self energy has the same dependence on the couplings and is 
already included. The lepton (or B-L) asymmetry thus produced then gets converted, through sphaleron 
processes, into the baryon asymmetry [9[, observed to be 0(10^^*^). In order to provide an acceptable 
framework for leptogenesis, and taking into account the presence of washout effects, the asymmetry we 
need is ei '-^ 0(10"^) or larger, for the same values of parameters that reproduce the observed small 
neutrino masses, Eq. ([¥]). In order to generate this asymmetry thermally one requires the Universe to 



reheat after inflation to a very high temperature 28|, |29|, |30|, |31 1 



Tr>2x 10^ GeV (14) 

Such large scale leads to an overproduction of cosmological gravitinos. In minimal supergravity models, 
with 7713/2 ~ 100 GcV to 10 TeV gravitinos are not stable, decaying during or after Big Bang Nucleosyn- 
thesis (BBN). Their rate of production can be so large that subsequent gravitino decays completely change 
the standard BBN scenario. Since the abundance of gravitinos is proportional to Tr such "gravitino cri- 



sis" can only be prevented by requiring a low enough reheat temperature Tr after inflation 10|, lllj . A 
recent detailed analysis gives a stringent upper bound 

Tr < 10^ GeV (15) 

when the gravitino has hadronic decay modes [l^l- Therefore, thermal leptogenesis seems difficult to 
reconcile with low energy supersymmetry if gravitino masses lie in the range suggested by the simplest 



minimal supergravity models. One possible way out is to have resonant leptogenesis jl3| as suggested 
in Ref. [14]. Another alternative considered in Ref. [ist requires going beyond the minimal seesaw by 
adding a small R-parity violating term in the superpotential. In the following we discuss quantitatively 
the alternative suggestion made in [16] in the context of the extended supersymmetric seesaw scheme. 



III. THE NEW EXTENDED SEESAW MODEL 

For definiteness we work in the context of the supersymmetric SO(IO) unified model considered in 
Ref.Q,[l3. The gauge symmetry and D-parity are broken by 45 and 210 multiplets at the unification 
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scale. The B-L symmetry breaks at lower scale thanks to expectation values of the x fields in the 16, 
instead of the more familiar left and right triplets present in 126 which lead to the standard seesaw 
mechanism. Gauge couplings unification can not fix the B-L breaking scale, which can be relatively low, 
as shown in Ref. [17]. The possibility of a low B-L breaking scale also fits with the observed neutrino 
masses. The relevant Yukawa couplings leading to neutrino masses are 



(16) 



Note that a direct Majorana mass term for the singlet fields Si is forbidden by an additional imposed 
U{1)g symmetry and the fact that the only singlet scalar present (cr) is odd under D-parity, while SiSj 
are even under D-parity. For the same reason, a cannot couple to S, but a bare mass is allowed Iq . 
As we will show leptogenesis allows this mass to be of the order of TeV. We also introduce a soft term 
breaking J7(1)g, which allows mixing between the scalar components of these fields ES'i, 



— £(3 = AjES'j, 

This will then give a 10 x 10 neutrino mass matrix, in the basis (z^^, S, v'[^, Si): 



(17) 



/ 





Yv 


FvL 


\ 





Ms 

















FvR 




^F^VL 


A 


F^VR 





/ 



(18) 



where, v = (0), vl — (xl) and vr = (xr) are the vevs for the fields 0, xl and xr respectively and A is 
the U{1)g breaking entry. Besides the three light neutrinos, this mass matrix will give two heavy states 
which are dominantly the right-handed neutrino v'^iL and the singlets Si, and a lighter state E. Notice 
that we will not impose any specific family symmetry. In Sec. IIV Bl we will, however, consider the case 
where the Yij are restricted to follow the Fritzsch texture for the quark sector. 



A. Light neutrinos 



Using the seesaw diagonalization prescription given in Ref. f26| we obtain the effective left-handed light 
neutrino mass matrix as 



1 

"My 



GG' - 



Y{FF-^f + {FF-')Y' 



vvl 
Vr 



(19) 



where G = Y{F~^)'^ A^. Since the D— parity breaking scale is much higher than the scale at which the 
left-right symmetry breaks, and this in turn is higher than the electroweak symmetry breaking scale, one 
has the "vev-seesaw" relation 



vl oc 



vrv 
Mx' 



(20) 



where Mx is determined by the SO(IO) breaking vevs. In the following we will take this proportionality 
constant equal to one and therefore write Eq. p9p as 



1 

Ms" 



GG' 



Y{FF-^f + {FF-')Y' 



V 

M^ 



(21) 
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Note that the neutrino masses are suppressed by the unification scale Mx instead of the B-L scale. As 
noted in 17[ we need two pairs of 16 16 in order to boost up the neutrino mass scale and bring in 
an independent flavour structure beyond that of the charged Dirac couplings, constrained by the quark 
sector phenomenology. 

Now we want to extract as much information from the low energy neutrino data as possible. The 
matrix mi, can be diagonalized as, 

U'^m^U = diag(TOi,^ , m^^ , m^^ ) (22) 
where we use the standard parameterization from Ref.jsl] and now adopted by the PDG 



U = 



( C12C13 S12C13 s-Y-ier"-^ 

-S12C23 - 0123235136"^ C12C23 - Si2S23Sl3e** S23C13 | (23) 
\ S12S23 - Cl2C23Sl3e*'' -C12S23 - Sl2C23Sl3e"' C23C13 



The values of U and di&%(jny^^my^^my^) are known to some degree from experiment and we want to 
turn this information into the parameters of the theory, defined by Eq. (|16p . Of course we have too little 
experimental information, three angles and two mass differences, to be able to reconstruct in full the 
three matrices Y, F, F and the vector A. 



B. Heavy neutrinos 



Let us now turn to the discussion of the heavy neutrinos. The following terms in the Lagrangian, 
= ^(2F,,iy',LSjVB, + MsJ:j: + 2A,J:S, 



IPs + VR S^ F' ty^L + A'^'ES* + AS' J: + MsSS 



lead to a 7 X 7 mass matrix in the basis N =(S, v^j^, Si): 

( Afs \ 



Mn = 



Fvr 
A F'^VR 



(24) 



(25) 



which determines the masses of the heavy neutrinos. Now in order to get the states with physical masses 
we diagonalize this mass matrix in three steps. We start by diagonalizing the matrix F . We define 



UlFUs VR = diag(Afi, M2, Af3) = Md 
where Un and Us are unitary matrices. Then we define new fields 



(26) 



= C/;^V^ S' = U^^S 



(27) 



so that in the basis N' ~ (S, v"^, S') we have 



-C = -N'^M'j^N' 



(28) 



7 



with 



( Ms A^f/s \ 
Md 
U^A Md 



(29) 



and 



Rns = 



f 1 \ 

Un 
Us 



(30) 



Now we rotate the fields ly"' and S' to obtain the almost degenerate states N± — -^{v"^ ± 5"). This is 
done by defining a new matrix 



R 



that obeys 



with N' = RN", N" = (S, N-,N+) and 



= R^M'^R : 



/l \ 
^2 \/2 

V° "71 71/ 



-£ = -N"^M'j(;N" 



(31) 



(32) 



Ms 



'Md 






Md 



(33) 



We now have to perform the final rotation in order to get the physical states for these heavy neutrinos. 
As we have seen in the discussion of the light neutrinos, the amount of G- violation must be small, therefore 
Eq. ((33)) can be approximately diagonalized using the techniques of Ref. 26| . We obtain 

1 



-C 



with N'" = (E', N'_,N'^), N" = UrhN'" where 



-N"'^M'j:;n"' 


V » 



(34) 






\ 


Md 








Md ^ 



(35) 



and 



Ub.h = 



-^M^'U^A 



V2 
1 



M^'U^A 



(36) 
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If we assume that the eigenvalues of F are hierarchical we get also a hierarchical spectrum for the heavy 
states, E', N'_,N'_^, namely, \M^\ <^ \Mi\ <^ IM2I <^ IM^]. In order to evaluate the asymmetry generated 
we must evaluate the couplings of the mass eigenstates. For this we notice that 



N = RnsRUrhN" 



(37) 



and a simple calculation gives 

/ 

RnsRUrh = 



1 





^/2 



%/2 



which allows us to write 



where we have used 



(38) 



(39) 



(40) 



With this we can rewrite the relevant part of the Lagrangian of Eq. in terms of the eigenstates 
(we drop the primes from now on). 



Cy = >s,: L,HY. + y±y N±, L.J H + h.c. 



(41) 



where 



(42) 



where an denotes the projection of the relevant light MSSM Higgs doublet h into the directions of the 
defining Higgs doublets living in H G 10^, and are the projections of the light MSSM-like Higgs 
doublet onto the defining Higgs doublets in the \Qfj and Mi is the mass of the N±i. 



C. Calculation of the Asymmetry 

We now discuss the issue of leptogenesis in this model. It can occur only after the local {B — L) C 
SO(IO) symmetry is broken. This will take place through the out-of-equilibrium decay of the singlet 
fermion E. The total width of E is given by (treating Ys as a column vector) 



~ -^FjtlsMs (43) 

OTT 



where Fe is given in Eq. 
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The asymmetry coming from the diagrams of Fig. [5] involves the sum over k which reduces to the sum 
over the hghtest pair of the (almost degenerate) N± = (y'^ ± S) states with masses ±Mi . A comparison 
of the diagrams in Fig. [5] and Fig. [T] gives the following dictionary 

Ajl — > Iei Ami Y^m , , ,s 

(44) 



We have then for the numerator 



V- .... 1 1 V- 



Yi^Yi„,A^B, (45) 



= -E 

where we have defined: Y±„i = A„i ± B„i. Comparing with Eq. (|42p we get 

Am = -^an {Y^ Um)^, , B,„ = --^a^ [FUs)^, (46) 
Putting everything together wc finally obtain 

3 M^ lH{YlFkUsa'^.)i{YlY^UNaH)i] 

SttMi yIy^ ^ ^ 

It is important to note that, in contrast to the asymmetry ei in the minimal seesaw discussed in the 
previous section, is not constrained by the light neutrino masses. This can essentially be understood 
by the fact that the neutrino mass, see eq. fTO]) and eq. ipO)) . is suppressed by the small ratio vl/vr, 
whereas in the calculation of e^, the small quantity A/w/? appears quadratically in the numerator and 
the denominator, and thus cancels, es can therefore be much larger than in the minimal seesaw case, 
independent of the light neutrino masses. Consequently, there is also no lower bound on A/s from the 
asymmetry parameter e^. 



IV. ANSATZES FOR THE COUPLING MATRICES 



We now estimate the resulting CP asymmetry needed for leptogenesis making use of the current values 
of the neutrino oscillsation parameters given in A simple ansatz is to assume that it comes just 
from the Dirac phase of the three-neutrino lepton mixing matrix assuming the unitary approximation. 
In this approximation the asymmetry is proportional to the unique CP invariant parameter that can be 




FIG. 2: Tree level and one loop diagrams for the decay of S that interferes to generate a lepton asymmetry of 
the universe. 
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probed in neutrino oscillations. The maximum value of the asymmetry that can be achieved is obtained 
by varying randomly all the other model parameters and the results are displayed at each point of the 
plane (sin0i3,sin i5). 



A. Generic case: non-hierarchical Yi 



We start by considering a non-hierarchical ansatz for the coupling matrices. In order to reduce the 
unknowns we assume that both Fij and Fij matrices are proportional to the standard SO(IO) Dirac 
Yukawa coupling matrix Yij, which is symmetric, namely 

With this choice Eq. (PT|) can be written as, 

= mt,'^) + m(,'') (49) 



where 



m 



= -ITT ^ ~3 ^ 10"' 1 (^0) 

Clearly, the projective nature of rrS^^ will not explain the neutrino data, therefore we need the contribution 



of rni^\ The atmospheric scale can easily be reproduced with Y ~ 0(1), /// — 0.1 and ~ 10 ^, for 
Afs = 1 TeV. 



With the assumptions of Eq. (|48p we can solve for the Yukawa matrix in terms of the experimentally 
observed neutrino oscillation parameters via 



1 



(51) 



which also involves other parameters of the model. 

With this ansatz and the choice / = 1, / = 0.1 we take the lepton mixing parameters in the range 
allowed by the experiment. At Scr the latest neutrino oscillation data give [l| 

sin^ 6*12 e [0.24, 0.40], sin^ 6*33 e [0.34, 0.68], sin^ 6*13 < 0.04 (52) 

In addition we take the other model parameters in the following ranges, 

A/s e [10^ 10^] GeV, A, e [10"^ 10^] GeV, vr, e [10^ 10^] GeV (53) 

With these values we can see the resulting w/j, Ms and As values that follow from our ansatz in Fig. [3J 
They show that it is possible to fit the neutrino data with this simple type of ansatz. The resulting CP 
asymmetry produced is given in Fig. |4l Here we have calculated the maximum value of the asymmetry 
that can be achieved at each point of the (sin ^13, sin (5) plane, varying randomly all the other parameters. 
Clearly the sizeable values of the asymmetry can be obtained, especially for large sin0i3 and sin (5, as 
expected, so that the necessary CP asymmetry needed for leptogenesis can be achieved. 

One sees that large values of the asymmetry, compatible with the experimental data, can be achieved 
in the full parameter space, even for very small sin0i3 and sin (5 values. However the ansatz is manifestly 
inconsistent with the successful Fritzsch texture for the quark masses. 
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0.00 0.05 0.10 0.15 0.20 

sine,3 

FIG. 4: Contour levels of the logarithm of the maximum value of the asymmetry obtained as a function of 
sin 013 and sin 5 for the non-hierarchical ansatz of section IV A. The other neutrino mixing parameters are taken 
the from the latest neutrino oscillation data given in 



B. 

It is natural to ask whether our SO(IO) model described by the Lagrangian of Eq. (fTB|) can provide 
thermal leptogenesis while reconciling the successful Fritzsch texture for the quark masses with the 
observed structure of lepton masses and mixings that follow from neutrino oscillation experiments. To 
this end we now assume that the Yij Yukawa's involved in neutrino mass generation are also restricted 
by the Fritzsch ansatz for the quark couplings, given in Eq. ((5|), with a, c complex and b,d real. Aware 
of the fact that the phases will be necessary in computing the final asymmetry, we will, for the moment 
take the ansatz parameters all real, as 



b= — , c=J — d= — 54 

These values imply a strong hierarchy among the Yukawa couplings. With this choice, let us now look at 
the neutrino mass matrix. It is clear from Eq. ||2T|) that this hierarchy must be corrected by a suitable 
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hierarchy in the F couphng matrix. In fact Eq. (PT|) can be "solved" for F as 

F = X^^F 



where X is obtained from 



YX' 



XY- = -^ 



u- 



YA 



2 

Me 



= Z 



(55) 



(56) 



where A = F^^A, F being, in general, an arbitrary non-symmetric matrix. Now Z is expressed as a 
combination of neutrino data and additional parameters, for which we take random values. 

We have performed a random study of this Fritzsch ansatz assuming as our ansatz that F ~ X~^F 
with F — i^"^ and taking random values of order one in various ways. We have found that in this case, 
for example for vr in the range [10^, 10*^] GeV, Ms of few TeV and As £ [IQ-^, 10*] GeV one indeed 
obtains a viable solution. One finds that some of the entries of F are small (of order 10^'^) in order to 
compensate for the corresponding smallness of Y in Eq. ([5|) . 

We have also explicitly calculated the value of the asymmetry es for this ansatz as explained above. 
The results are shown in Fig. [5] for the case of the hierarchical ansatz with the current neutrino oscillation 
parameters from ^I'j and the Dirac phase is (5 e [0, tt]. 
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FIG. 5: Contour levels of the maximum asymmetry obtained as a function of sin 6*13 and sin 5 using the Fritzsch 
texture with F — X^^F, F=I (left panel) or F symmetric (right panel). 

One sees that, in this case, values of the asymmetry of order 10^^-10^^ can only be obtained for large 
values of sin ^13 and sin (5. Therefore, this scenario can be potentially probed by the future neutrino 
oscillation measurements. 

Larger values of the CP asymmetry compatible with current neutrino oscillation measurements can 
be found for other ansatze of the current model using the Fritzsch texture. For instance, if we consider 
F = fY with F symmetric, or F diagonal and F — XF, one can obtain very large values of the asymmetry 
even for very small values of sin ^13 and sin J, as illustrated in Fig. [6l 



V. ASYMMETRY WASHOUT, SINGLET PRODUCTION AND THE SPHALERON 

CONSTRAINT 

Generating a large enough asymmetry in the decay of the lightest singlet is a necessary, but not 
a sufficient condition for successful leptogenesis. Additional conditions are required, before one can 
conclude that any given model can generate the required baryon asymmetry. 
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The first condition to be satisfied is the out-of-equilibrium decay of the heavy singlet, which is nothing 
but the basic Sakharov condition [32]. This requires that the decay rate is smaller than the expansion 
rate of the universe, i.e. Ts ^ THubbic at the decay epoch. For large Fx;, say, of order of the expansion 
rate, part of the asymmetry produced in the decays will be washed out by inverse scattering processes 
violating lepton number. This constraint will thus put an upper limit on the couplings of S to the leptons 
(and Higgses) for any given mass M^- Second, we need to produce a sufficient number of singlets in the 
early universe. Singlets could be either produced by their couplings to the thermal bath or through the 
decay of some heavier particle, which was present in the universe at earlier times. Obviously, a sufficiently 
strong direct thermal production requires a minimum value for the couplings of S, whereas the second 
option does not. The third constraint comes from the fact that the S decays have to take place at a time 
early enough that the SM sphalerons are still in equilibrium, otherwise we will produce a lepton number, 
but no non-zero baryon number. This constraint puts a lower limit on Fs, F^ > Fgphaicron, independent 
of the production mechanism of E. 

To accurately calculate the first two of the above conditions, in principle, one needs to set up a 
network of Boltzmann equations, which in general can only be solved numerically [s^. However, under 
certain simplifying assumptions, one can derive approximate analytical solutions which reproduce the full, 
numerical calculations quite well. Several analytic al ap proximations have been proposed in the literature, 

|. In our earlier paper [l6| we have used the following approximate form for 



33 



see for example 



the washout factor: 



34 



k'{z) = - 



102 



(57) 



where z = f—^^-^: with Ffjubbic being the expansion rate of the universe 33 1. Ref. 35 1 has numerically 
solved the Boltzmann system for the case of the simpler type-I seesaw. In this case the decay width of the 
lightest right-handed neutrino is proportional to F ~ {Y^Y)ii, and the authors of [35[ define the "effective 
mass" parameter mi = (F^y)iif ^/Mi, where Mi is the mass of the lightest right-handed neutrino and 
V the SM vev. Thermal equilibrium for the right-handed neutrino is reached for an "equilibrium mass" 
m^: = 1.08 • 10^'^ eV, for which by definition z — 1. The full numerical calculation is then very well 
approximated by 



k(z) ~ 0.24(a;_ 



(58) 
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with 



and the numerical values 



x± 



Vm±/ 



(59) 



B.3 X 10"'' eV, 



3.5xlO-''eV, a = 0.1. 



(60) 



(61) 



We can make use of the results of with the straightforward replacement, 

x± = ( 1 

Our more complicated setup of Boltzmann equations then is solved approximately by the fit eq. (|58p . 
Note, that this implicitly assumes (40|] that the heavier singlets are sufficiently decoupled so as to not 
contribute significantly to the washout. 

As seen in Fig. ([7]) eq. ((58|) and eq. ([57|) lead to very similar results for z > 1, The two forms 
differ, however, for z < 1. This can be traced to the fact that eq. (|58p also accounts for the suppressed 
production in the weak coupling regime, while eq. (|57p does not. For definiteness, in the plots shown 
below, we have used eq. (fSH]) . 




FIG. 7: Calculated range for r]B as a function of z, for the different fits to the washout factor. The dashed (full) 
lines show the approximate analytical solutions to the Boltzmann equations, see eqs (I57|l and (|58|) . Results are 
shown for 2 typical values of e = 10~^, 10~®. The horizontal black band is the currently allowed r)B range [stI ]. 



We have also used the constraint imposed by the condition that the decay of the singlets must happen 
while the sphalerons are still active. We estimate the sphaleron time using 

y2 



Tsphalcron — ^ YlnhhlcW = EW PT — \/ 47r-^5*/45(— — y\T=EWPT-, (62) 

with EWPT denoting the energy at which the electro- weak phase transition occurs, the effective 
number of degrees of freedom and Mpi the reduced Planck mass. We cut all points with Fg < Fgphaieron- 
This is, of course, a rough approximation to the real, dynamical situation. However, we believe it to be 
sufficiently accurate for deriving order of magnitude constraints on the model parameters. Finally, when 
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converting the produced lepton asymmetry to the baryon asymmetry, we have to take into account an 
efficiency factor for the sphalerons. This factor has been calculated in [36[ to be 



VB 



H 



22nF + ISuH 



(63) 



where np {nn) is the number of families (Higgses). Numerically this factor is ^ 1/3. 

Fig. ([8]) shows as an example the resulting r]B as a function of A for a numerical scan using the ansatz 
shown in fig. ([5]), to the left. The light/dark (green/red) area is the calculated range for r/s without/ with 
the sphaleron constraint. One sees that the different constraints discussed above conspire to choose a 
rather well-defined allowed range for the parameter A in this case. A large enough baryon asymmetry 
can be obtained roughly for A = [1, 10^] GeV, for the range of the other parameters as given in eq. ((53|) . 
It is amusing to note, that the requirement of producing a sufficient number of singlets and the sphaleron 
constraint lead to rather similar lower cuts on A 



41| 




10 10 10 10 10 10 10 10 10 



A [GeV] 



FIG. 8: Calculated range for riB as a function of A, for the ansatz defined in fig. ([5]), to the left. The horizontal 
band is the allowed range of rjB- The full shaded region corresponds to the allowed range without the sphaleron 
constraint, while the superimposed darker/red area indicates the allowed range of parameters once the sphaleron 
constraint is included. For a discussion see text. 



We have repeated this exercise for all the different ansatze discussed above. The resulting allowed 
ranges for the parameters A and are shown in fig. ([9]). As shown in this figure, the random non- 
hierarchical ansatz and the hierarchical ansatz shown in the right panel of Fig. [6] lead to the widest 
allowed ranges for A and Ms- 

Before we close this section let us briefly illustrate possible values for the relevant parameters. These are 
A, Afs, the masses Mi, M2, M3 of the three iso-singlet neutrinos and the resulting baryon asymmetry 
rjB- Some examples are given in table |T] for each of the scenarios we have discussed. As one can see the 
acceptable baryon asymmetry may arise for rather low scales, in contrast to the (simpler) seesaw type-I 
schemes. Note, however that the values given in the table are by no means unique, and are also not 
meant to be "representative" of the classes. The generation of the baryon asymmetry is far easier in 
this context than in the traditionl seesaw model. The presence of the new singlets allows us to have, in 
addition, acceptable textures for the quark and lepton mixing angles. 
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Non-hierarchical Hierar. 6b 



FIG. 9: Calculated ranges for A and A/e for the different texture ansatze discussed in the text. This ranges have 
been obtained with the largest possible scan for the parameters, for example, Me in the full interval [10^, 10^^] 
GeV. 



A[GeV] 


Me [GeV] Mi [GeV] A/2 [GeV] M3[GeV] 77s [10"'°] 


Non-hierarchical 


8.3x10-2 


1.2x10^ 


4.1x10* 5.0x10* 


2.6x10^ 


6.3 


9.8 


2.3x10* 


9.3x10^ 1.2x10^ 


6.1x10^ 


5.0 


Hierarchical 5a 


3.5x10^ 


1.2x10^ 


3.6x10* 6.2x10* 


6.7x10® 


4.9 


2.9x10^ 


1.1x10^ 


4.6x10* 6.8x10* 


5.7x10*^ 


5.8 


Hierarchical 5b 


1.8x10^ 


6.6x10^ 


2.1x10* 6.4x10* 


3.8x10® 


5.2 


8.3x10^ 


2.5x10^ 


6.1x10^ 6.2x10* 


9.8x10® 


4.8 


Hierarchical 6a 


2.9x10^ 


2.9x10^ 


1.6x10^ 3.6x10^ 


6.0x10" 


5.1 


3.3x10^ 


8.3x10^ 


2.5x10^5 1.3x10" 


1.4x10" 


5.8 


Hierarchical 6b 


1.0x10^ 


1.8x10^ 


1.9x10^ 1.1x10'^ 


2.0x10® 


5.4 


3.4x10^ 


2.8x10^ 


1.5x10'^ 1.5x10'^ 


2.5x10® 


5.9 



TABLE I: Some example spectra for the different ansatze considered. 
VI. CONCLUSIONS 

We have reviewed the argument that in minimal type-I SO(IO) seesaw one can not easily reconcile 
the thermal leptogenesis scenario with the successful Fritzsch texture for the quarks and an acceptable 
pattern of lepton masses and mixings that follow from neutrino oscillation experiments. This is due to the 
fact that the large seesaw scale needed to account for small neutrino masses leads to an overproduction 
of cosmological gravitinos, which destroys the standard predictions of Big Bang Nucleosynthesis (BBN). 
Barring the very special case of resonant leptogenesis, one must go beyond the minimal type-I seesaw 
mechanism. 

In this paper we have studied in some detail an extended seesaw scenario as a natural way to overcome 
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this limitation. The proposed extension has the added virtue of providing a natural setting for reconciling 
the observed structure of neutrino mixing angles with the strong hierarchy among quark masses and the 
smallness of the quark mixings. While this can always be accommodated in a "generic" unified theory 
with arbitrary multiplet content, it becomes a real challenge for unified predictive models of fiavour. 

We have provided a quantitative study of fermion masses and (an approximative calculation of) lep- 
togenesis in the context of supersymmetric SO(IO) Unification. Our approach was phenomenological in 
that we have not assumed a specific flavour symmetry. We have shown how thermal leptogenesis can 
occur at relatively low scale through the decay of a new singlet, thereby avoiding the gravitino crisis. 
Washout of the asymmetry is effectively suppressed by the absence of direct couplings of the singlet to 
leptons. For illustration we have shown how one can accommodate current oscillation neutrino data and 
the required value for the asymmetry for successful leptogenesis even if the only source of CP violation 
is the Dirac phase S in the low energy neutrino mixing matrix. Finally, we note that we have not taken 
into account flavour effects 33, |39| in the calculation of the asymmetry. However, we believe that the 



absence of a lower bound on Ms in our model is independent on whether flavour effects are taken into 
account or not. 

Using the Fritzsch texture we have found that some ansatze lead to acceptable values of the asymmetry 
of order 10~^-10^^ only for large values of sin 6*13 and sin (5. Therefore, these scenarios can be potentially 
probed by the future neutrino oscillation measurements, as the required value of the CP invariant is nearly 
maximal. In contrast, we have also presented alternative Fritzsch-type ansatze leading to substantially 
larger values of the CP asymmetry, even for very small values of sin ^13 and sin 6. We have also discussed, 
how to approximately treat the conversion of the decay asymmetry to the baryon asymmetry, without 
resorting to a full numerical solution of the Boltzmann equations. To this end, we made use of some 
approximation formulas derived previously for seesaw type-I and discussed how they can be adapted to 
cover also our more complicated case. 

In summary, our extended seesaw scenario provides a way of reconciling the lepton and quark mixing 
angles with thernal leptogenesis in a unified scenario. While this by itself does not constitute a complete 
theory of fermion masses and leptogenesis, at least it provides a useful first step towards an ultimate 
unified theory incorporating flavour. 
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